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ABSTRACT 

We calculate the tidal response of helium and carbon/oxygen (C/O) white dwarf (WD) binaries inspiraling 
due to gravitational wave emission. We show that resonance locks, previously considered in binaries with 
an early-type star, occur universally in WD binaries. In a resonance lock, the orbital and spin frequencies 
evolve in lockstep, so that the tidal forcing frequency is approximately constant and a particular normal mode 
remains resonant, producing efficient tidal dissipation and nearly synchronous rotation. Resonance locks in 
WDs can occur not only with global standing modes, but even when damping is so efficient that the resonant 
tidal response becomes a traveling wave. We derive simple analytic formulas for the tidal quality factor Q t and 
tidal heating rate during a g-mode resonance lock, and verify our results numerically. We find that Q t ~ 10 7 for 
orbital periods < 1 - 2 hr in C/O WDs, and Q t ~ 10 9 for P or b < 3 - 10 hr in helium WDs. Typically tidal heating 
occurs sufficiently close to the surface that the energy should be observable as surface emission. Moreover, near 
an orbital period of ~ 10 min, the tidal heating rate reaches ~ 10~ 2 L Q , rivaling the luminosities of our fiducial 
WD models. Recent observations of the 13-minute double- WD binary J0651 are roughly consistent with our 
theoretical predictions. Tides naturally tend to generate differential rotation; however, we show that the fossil 
magnetic field strength of a typical WD can maintain solid-body rotation down to at least P or b ~ 10 min even 
in the presence of a tidal torque concentrated near the WD surface. 
Subject headings: binaries, white dwarfs, waves, instabilities 



1. INTRODUCTION 

In this work, we consider the effect of tides in detached 
white dwarf binaries inspiraling due to energy and angular 
momentum loss by gravitational waves. Our analysis is mo- 
tivated by several important questions. For example, to what 
degree should short-period white dwarf binaries exhibit syn- 
chronized rotational and orbital motion? Should white dwarfs 
in close binaries be systematically hotter than their isolated 
counterparts, as a result of tidal dissipation? What is the ther- 
mal state of white dwarfs prior to the onset of mass transfer? 

Several past studies have applied linear perturbation the- 
ory to the problem co nsidered in this work. Campbell ( 1984 ) 
and |Iben et aT] ( |1998| l applied the theory of the equilibrium 
tide to white dwarf binaries , using parameterized viscosities 
to estimate the tidal torque. Wille ms et al.| p010) considered 
turbulent convective damping acting onme equilibrium tide, 
as originally considered by Zahn ( 1977 1 for late -type stars, 
and showed that this effect is not able to synchronize a white 
dwa rf binary within its gravitatio nal wave i nspir al time. 

|Rathore et al~| ( |2005] > and puller & Lai] ( pOlT) moved be- 
yond the large-scale, nonresonant equilibrium tide, and con- 
sidered the tidal excitation of standing g-modes during inspi- 
ral, analyzing the behavior of wave amplitudes as a system 
sweeps through resonances. However, neither study allowed 
the white dwarfs spin rate to evolve, an assumption that elim- 
inated the physics highlighted in this work. 

In this paper we also focus on tidally excited g-modes in 
white dwarf binaries; one of our goals is to assess whether 
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the resonantly excited "dynamical tide" represents a traveling 
or standing wave. This amounts to whether a tidally gener- 
ated wave is able to reflect at its inner and outer radial turn- 
ing points. If reflection cannot occur, then a damping time of 
order the group travel time across the mode propagation cav- 
ity results; if reflection does occur, then the wave amplitude 
can build up significantly during close resonances between g- 
mode frequencies and the tidal forcing frequency. 

This question has b een addre ssed before in the context of 
main-sequence stars. Zahn|( [T975] l employed a traveling wave 
description of the dynamical tide in the context of early -type 
stars with radiative envelopes, assuming waves would be ab- 



sorbed near the surface by rapid radiative diffusion. Goldre 



ich & Nicholson (1989a) enhanced this argument, showing 
that dynamical tides first cause tidal synchronization in such a 
star's outer regions, leading to the development of critical lay- 
ers and even stronger radiative damping. However, they did 
not assess whether angular momentum redistribution could 
enforce solid-body rotation and thereby e limin ate critical lay- 
ers; we address this imp ortant point in § 6.2 In the absence 
of critical layers, Witte & Savonije ( 1999) introduced the phe- 
nomenon of resonance locks (§ |4J, which rely on the large 
wave amplitudes produced during eigenmode resonances. We 
will show that such resonance locks occur ubiquitously in 
clos e white dwarf binarie s. 

Goodman & Dickson (1998 ) considered the case of late- 
type stars with convective envelopes, and showed that tidally 
generated waves excited at the edge of the conv ection zone 
steepe n and break near the cores of such stars. |Fuller & Lai| 
(2012), in their study of the tidal evolution of white dwarf 



binaries, found that the dynamical tide in a carbon/oxygen 
white dwarf instead breaks near the outer turning point. As 
such, they invoked an outgoing-wave boundary condition in 
their analysis. We find that their assumption may not be gen- 
erally applicabl e du e to an overestimate of the degree of wave 
breaking; see § 6.1 As a result the dynamical tide may repre- 
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sent a standing wave for a substantial portion of a white dwarf 
binary's inspiral epoch. 

This paper is organized as follows. In §[2] we provide per- 
tinent background information on white dwarfs and tidal ef- 
fects that our subsequent results rely on. In § [3] we give a 
broad overview of the results we derive in more detail in §§|4] 
-|8] In §|4]we consider the case of resonance locks created by 
standing waves. We analyze the resulting tidal efficiency and 
energetics in §[5] In §[6] we analyze whether sta ndin g waves 
i occur, considering wave br 
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are able to occ ur, co nsidering wave breaking in § |6.1| and crit- 
ical layers in § |6.2| In §[7] we turn our attention to trave ling 
waves, discussing wave excitation and interference in § |7.1| 
and showi ng t hat traveling waves can also create resonance 
locks in § 7.2 In §[8] we then employ numerical simulations 
to combine our standing and traveling wave results. In §[9] we 
compare our results to observational constraints and discuss 
physical effects that need to be considered in future work. We 
then conclude in § 10 with a summary of our salient results. 



2. BACKGROUND 

A short-period compact object binary efficiently emits grav- 
itational waves that carry off energy and angular momen- 
tum. This process causes its orbit to circularize; as such, 
we will restrict our attention t o circular white dwarf binaries 
in this work (however, see |Thompson||201 1) . Gravitational 
waves also cause such a binary's orbit to decay according to 
ft = f2/i gw , where the characteristic gravitational wave inspi- 
ral timescale for a circular orbit is (Pet ers|1964[ ) 



5 (l+M'/M) 
96 M'/M 
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-ft 
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(1) 



Here u> \ = GM/R 3 is the dynamical frequency of the primary, 
M' is the mass of the companion, 01 = GM/Rc is the relativ- 
ity parameter of the primary, and £1 is the orbital frequency. 
The time a binary will take until it begins to transfer mass 
is given by f merge = 3f gw /8. For a 0.6M Q white dwarf with an 
equal-mass companion, orbital periods of less than ~ 530 min 
imply the binary will begin mass transfer within 10 Gyr; this 
restriction reduces to P or b < 270 min for a 0.2M Q white dwarf 
with an equal-mass companion. 

During the process of inspiral, the tidal force acting on each 
element of the binary steadily grows. The tidal response on a 
star is typically divided conceptually into two components: 
the equilibrium tide and the dynamical tide. The equilibrium 
tide represents the l arge-scale distortion of a star by a com- 
panion's tidal force (Zahn 1977); it is often theoretically mod- 
eled as the filling of an equipotential surface, but can also be 
treated as the collective nonresonant response of all of a star's 
eigenmodes. The two viewpoints are equivalent, as in both the 
tidal forcing frequency is set to zero. Except near very strong 
resonances, the equilibrium tide contains the great majority 
of the tidal energy. Nonetheless, whether it produces a strong 
torque is also influenced by the degree to which it lags behind 
the tidal potential, which is determined by ho w strongly the 
equilibrium tide is damped. For white dwarfs, Willems et al. 
( |2010| ) showed that turbulent convection acting on the equi- 
libriu m tide does not cause significant synchronization (Ap- 
pendix |BT2]i. 

The dynamical tide, on the other han d, correspon ds to the 
tidal excitation of internal stellar waves (Za hn|1975| l. In par- 
ticular, given that tidal forcing periods are much longer than 
the stellar dynamical timescale, buoyancy-supported gravity 
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FIG . 1 . — Propagation diagrams for several of our white dwarf models listed 
in Table^ The top two panels are our 0.2Mq He 10 and He 5 helium white 
dwarfs, with T c f[ = 9,900 and 5,100 K respectively; the bottom two panels are 
our 0.6M Q C012 and C06 carbon/oxygen white dwarfs, with T efl = 12,000 
and 5,500 K respectively. Each plot shows the Brunt- Vaisala frequency N 
(green line; dashed indicates N 2 < 0), the quadrupolar Lamb frequency S 
(thick blue line), and inverse local thermal time l/f t h = gF /pc p T (red dot- 
dashed line). A g-mode is able to propagate where its frequency is less than 
both N and S. In the bottom panel showing our CO 6 model, the shaded region 
at high pressure indicates the plasma interaction parameter T > 2 20, im plying 
crystallization occurs (which is not included in our model); see §|9.3| 



waves or g-modes are predominantly excited (although rota- 
tionally supported modes become important when the rot ation 
and tidal forcing frequencies become comparable; see § 9.2 1. 
Propagation of gravity waves is primarily determined by the 
Brunt-Vaisala frequency /V, given by 



1/1 dlnp dlnp 



(2) 



dr dr 

where Y\ is the adiabatic index. A g-mode is able to propagate 
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TABLE 1 
White dwarf models 



ID 


Teff (K) 


L/Lq 




R/Rq 


f* (sec) 


fthlrcb (F) 


tool (Gyr) 


M com /M 


Tcore 


/W10- 2 


h/ MR 2 


HelO 


9,900 


1.3 x 10" 


-2 


0.038 


26. 


7.7 x 10 -8 


0.36 


3 X 10~ 14 


1.2 


0.34 


0.085 


He7 


7,000 


1.9 x 10" 


-3 


0.029 


18. 


5.1 x 10 1 


0.98 


1 x 10~ 7 


2.7 


0.38 


0.11 


He5 


5,100 


3.9 x 10" 


-4 


0.025 


14. 


1.3 x 10 6 


2.5 


2x 1Q- 4 


5.1 


0.41 


0.14 


C012 


12,000 


3.0 x 10" 


-3 


0.013 


3.1 


6.5 x 10 -9 


0.59 


1 x 10 -16 


71. 


0.99 


0.16 


C06 


5,500 


1.3 x 10" 


-4 


0.013 


2.9 


7.3 X 10 2 


4.5 


3 x 10~ 8 


260. 


1.0 


0.18 



Masses of the helium models (top three) are 0.2Afo ; masses of the carb on/oxygen models (bottom two) are 0.6Mq. Helium models each have a hydrogen layer 
of mass 0.0033M, and were generated with MESA (Paxton et al. 201 1 1; carbon/oxygen models each have a helium layer of mass 0.017M and a hydrogen layer 
of mass 0.0017M. Further details on white dwarf models are given in Appendix |C| Figure[T]provides propagation diagrams for several models. The dynamical 
time is t\ = S 3 /GM; the thermal time at the radiative-convective boundary (RCB) is ? t hlrcb = P c pT/gF | rc |,, where we take 2ir/N = 100 min to define the RCB; 
the white dwarf cooling time is f coo i = E^/L, where Efa = f c p TdM approximates the total thermal energy; M con v is the mass in the outer convection zone, 
which increases in size by many orders of magnitude as a white dwarf cools (see Figure[TJ; the plasma interaction parameter at the center of the white dwarf is 
r col - e = Z 2 e 2 /kTdi\ r= o, where Ze is the mean ion charge, d x is the ion spacing, and the value of T corresponding to the onset of crystallization is discussed in 
§|9.3| the relativity parameter is /3 2 = GM/Rc 1 ; and /* is the white dwarf moment of inertia. 



where its frequency lies below both N as well as the Lamb 
frequency Sj = 1(1 + l)c 2 jr 2 . Plots of both N and Si for several 
helium and carbon/oxygen white dwarf models are provided 
in Figure [T] 

Degeneracy pressure satisfies p oc p r ' ; substituting this into 
equation |2]) yields N = 0. Thus the Brunt-Vaisala frequency 
becomes very small in the white dwarf core where degener- 
acy pressure dominates gas pressure, scaling as A^ 2 oc kT /Ep, 
where E F is the Fermi energy. Moreover, white dwarfs also of- 
ten possess outer convection zones with N 2 < 0. As a result, a 
typical tidally excited g-mode in a white dwarf possesses both 
an inner turning point near the core as well as an outer turning 
point near the radiative-convective boundary. 

Lastly, temporarily ignoring degeneracy pressure and as- 
suming an ideal gas equation o f state, equation |2) can be ex- 
pressed as ( |Hansen et al. 12004] ) 



N 2 = ad/ ^ " (V ad -V)-g 



dr ' 



(3) 



where // is the mean molecular weight. From this expression 
we can see that the Brunt-Vaisala frequency becomes larger 
in composition gradient zones, where [i decreases with radius. 
This can be seen in Figure [T| where a "bump" in N occurs in 
helium models due to the helium to hydrogen transition; two 
bumps are present in carbon/oxygen models, resulting from 
carbon/oxygen — > helium and helium — > hydrogen. 

3. DYNAMICAL TIDE REGIMES IN WHITE DWARFS 

Here we will give a general overview of the results covered 
in §§F] _ [l]t , y enumerating four essential regimes of the dy- 
namical tide in white dwarfs, which comprises the wavelike 
tidal response. The two basic distinctions made by our four 
regimes are a) whether tidally excited gravity waves can re- 
flect and become large-amplitude standing modes, or whether 
they instead represent traveling waves; and b) whether or not a 
resonance lock ca n be created. Resonance locks are described 
in detail in §§|4|& [772{ they occur when the tidal torque causes 
the tidal forcing frequency a = 2{fl-H, spm ) to remain constant 
even as the orbit shrinks. 

First, the two regimes where standing waves exclusively oc- 



Sl) In this regime the dynamical tide represents a purely 
standing wave, but with a resulting torque that is in- 
sufficient to effect a resonance lock even during a 



perfect resonance. This occurs for long orbital peri- 
ods or small companion masses. As such, the sys- 
tem quickly sweeps through resonances, and the time- 
averaged torque is dominated by its value away from 
resonances. This nonresonant torque is proportional to 
the eigenmode damping rate, and is thus very small for 
white dwarfs, due to their long thermal times. As a con- 
sequence the average tidal quality factor is very large, 
tidal heating is negligible, and the spin rate remains es- 
sentially constant. 

S2) Here the dynamical tide is again a standing wave, but 
with eigenmode amplitudes large enough to create res- 
onance locks. This regime is addressed in detail in §|4j 
we estimate the orbital period corresponding to its on- 
set in equation ( fT0| ). During a resonance lock, tides 
become efficient: due to strong tidal torques, the spin 
frequency changes at the same rate as the orbit decays 
due to gravitation wave emission. Definite predictions 
result for the tidal energy deposition rate (equation 20 1 
and tidal quality factor (equation 17 1. 



Next, the regimes strongly influenced by traveling waves 
are: 

Tl) In this regime, the off-resonance dynamical tide is still a 
standing wave, but near resonances the wave amplitude 
becomes so large that reflection n ear th e surface can- 
not occur due to wave breaking (§ 6.1 1. Furthermore, 
the traveling wave torque is too weak to cause a reso- 
nance lock. Since the typical torque experienced by the 
white dwarf is once again the off-resonance standing 
wave value, the synchronization and tidal heating sce- 
nario is very similar to regime (SI) — in other words, 
tides are ineffective. 

T2) Just as with regime (Tl), the standing wave torque is 
capped at resonances, becoming a traveling wave; how- 
ever in this regime the traveling wave torque itself is 
strong enough to create a resonance lock, as addressed 
in § [7| W e estimate the onset of this regime in equa- 
tion ( p3[ ). The predictions for the tidal energy deposi- 
tion rate and quality factor are the same as in (S2). 

Although we consider only these four regimes in this work, 
at shorter orbital periods and nearly synchronous rotation, 
physical effects such as Coriolis modification of stellar eigen- 
modes and nonlinear tidal exc itatio n mechanisms are likely to 
become very important; see § |9.2| 
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The archetypal scenario is as follows. A white dwarf bi- 
nary with a sufficiently long orbital period begins in regime 
(SI). Eventually, as the orbit shrinks due to gravitational ra- 
diation and the tidal force correspondingly increases in mag- 
nitude, the dynamical tide becomes strong enough that a res- 
onance lock takes effect and regime (S2) is reached. How- 
ever, as inspiral accelerates, the torque necessary to maintain 
the resonance lock becomes larger, requiring larger wave am- 
plitudes. When the amplitude becomes too great, the wave 
begins to break near the outer turning point, and the system 
enters regime (Tl). Finally, when the traveling wave torque 
becomes large enough to create a resonance lock, it enters 
regime (T2). In §[8] we verify this picture numerically. 

4. RESONANCE LOCKS BY STANDING WAVES 

We assume in this section that the dynamical tide is a su- 
perposition of standing waves and proceed to predict the tidal 
evolution of a white dwarf binary. We discuss the applicabil- 
ity of the standing wave limit in § [6] 

Assuming a circular orbit and alignment of spin and orbital 
angular momenta, the secular tidal torque on a star can be ex- 
pressed as a sum over quadrupolar (/ = 2) eig enmo des indexed 
by their number of radial nodes n (Appendix B.2 1: 



t = &mE*e 



(oj 2 -a 2 ) 2 +4a 2 j 2 



(4) 



Here a = m{Vt - O sp j n ) is the tidal driving frequency in the 
corotating frame, 51 is the orbital frequency, O sp ; n is the solid- 
body rotation rate, m = 2, = GM 2 /R is the white dwarf 
energy scale, e = (M' /M)(R/a) 3 is the tidal factor, M' is the 
companion mass, a is the orbital separation, W 2 = 37r/10, and 
7„ is an eigenmode damping rate ( Appendix |B.3| >. Our e igen - 
function normalization convention is given in equation ( |B5[ ); 
physical quantities such as the torque are of course indepen- 
dent of the choice of normalization. 

The linear overlap integral Q„ appearing in equation Q 
represents the spatial coupling strength of an eigenmode to 
the tidal potential, and is normalization dependent. Since the 
tidal potential spatially varies only gradually, Q„ is large for 
low-order modes, and becomes much smaller for high-order, 
short- wavelength modes. We describe various methods of 



computing Q„ in Appendix B.4 



The factor in brackets in equation (|4j) describes the tem- 
poral coupling of an eigenmode to the tidal potential, and 
becomes very large during resonances, when the tidal driv- 
ing frequency becomes close to a stellar eigenfrequency. The 
nonresonant limit of this factor, which corresponds to the 
equilibrium tide, increases with stronger damping. Paradoxi- 
cally, however, the torque during a resonance is inversely pro- 
portional to the damping rate, since damping limits the maxi- 
mum energy a resonant mode attains. 

We note that by invoking steady-state solutions to the mode 
amplitude equations, as we have done here, we fail to account 
for the energy and angular momentum transfer required to 
bring a mode's amplitude up to the steady-state value. Ad- 
ditionally, the steady-state solution itself may fail to model 
the behavior of mode amplit udes very close to resonances 
correctly; we address this in § |6.3| Correctly accounting for 
these two considerations would involve simultaneously solv- 
ing both the mode amplitude equations for all relevant modes 
as well as the orbital evolution equations, a task we leave to 
future study. 



Continuing, we focus on resonant tidal effects, and consider 
the case of a particular eigenmode with a frequency close to 
the tidal driving frequency, i.e. lj„ ss a. We can then make 
this substitution everywhere in equation Q other than in the 
detuning frequency Soj„ = uj„ - a to find 



t « 2mE„e l W 2 Ql 



(5) 



having dropped nonresonant terms. One might expect that 
since the strongest torques are achieved very near resonance, 
where 6lj„ ~ 0, a system should evolve quickly through reso- 
nances, and they should have little effect on the long-term or- 
bital and spin evolution. This is often accounted for by using a 
"harmonic mean" of the torque to produce a synchronization 
time, i.e., f svnc = f I^dfl/r, w here is the moment of inertia 
( [Goodman & Dickson|1998] l. 

Under particular circumstances, however, it is possible to 
achieve a resonance lock, where an eigenmode remains in a 
highly resonant state for an extended period of time, as origi- 
nally proposed by |Witte & Savonije"| ( |1999| ). Very near a res- 
onance, the torque depends very strongly on the detuning fre- 
quency 5uj„, and very weakly on the orbital period by itself; 
as a result, the essential criterion that must be satisfied for a 
resonance lock to occur is that the detuning frequency must 
remain constant: 



= 5u)„ ■■ 



,((.- 



C "lO 



spin 



(6) 



where dui„ / <9il sp i n = -mC n accounts for rotational modifica- 
tion of the stellar eigenmodes , C n ~ 1/6 for h igh-order / = 2 
g-modes and slow rotation ( |Unno et al.||1989| l, and we have 
assumed the white dwarf rotates as a s olid body. (We justify 
the solid-body rotation assumption in § |6.2| ) 

For simplicity, we will henceforth ignore rotational modifi- 
cation of the stellar eigenfrequencies, so that C„ — > 0. This 
limits the quantitative applicability of our results to where 
^-^spin ^ ^spin- We also neglect progressive white dwarf 
cooling, which decreases the Brunt-Vaisala frequency and 
consequently lowers eigenmode frequencies; this is valid so 
long as the cooling time f coo i, which is on the order of ~ Gyr 
for the models listed in TablefT] is much longer than the grav- 
itational wave decay time f g „(equation[T]). Subject to these 

simplifications, equation |6]) then reduces to Cl = fl sp - m , i.e., 
that the orbital and spin frequencies evolve at the same rate. 
Since the orbital frequency increases due to the emission of 
gravitational waves, and the spin frequency increases due to 
tidal synchronization, this phenomenon is plausible at first 
glance. We now work out the mathematical details. 
The evolution of f2 sp i n and £1 proceed as 



a 



spin / ^spin 

n/n 



l/? gw + (3/2)(£ tide /|£ orb |) 



(7) 



where E^ is the secular tidal energy transfer ratq^j and the 
gravitational wave inspiral time f gw is given in equation (TT]|. 
Here we have failed to account for tidal effects in the com- 
panion, which would provide an extra contribution to Ex x & t \ 
see below. 

Using equations (|7]i and ( |B8| ), equation |6]l becomes 



n 



l 



3 1 

2 fia 2 



(8) 



5 Our convention is that r > or £" t ide > implies that orbital angular 
momentum or energy is being transfered to the white dwarf(s). 
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FIG. 2. — Plot schematically illustrating the dynamics of a resonance lock. 
The abscissa is the orbital frequency Q, which increases due to gravitational 
wave radiation, and the ordinate is the Doppler-shifted I = m = 2 tidal driv- 
ing frequency a — 2(f2 — n spm ). The arrows depict the vector field describing 
the orbital evolution equations (equation [7). The eigenfrequency of the in- 
cluded mode is uj, which is flanked by stable and unstable fixed points of 
the evolution equations: uj- and uj+, respectively. The dashed red horizontal 
lines show these three frequencies, while the blue curves are example sys- 
tem trajectories. The stable point ui- corresponds to a resonance lock, while 
the unstable fixed point u + corresponds to the upper boundary of the stable 
fixed point's basin of attraction (shaded region). In producing this plot, the 
correct functional form of the equation of motion has been used, except with 
f gw and e taken as constant, and with artificially chosen values of the various 
germane parameters. In particular, for realistic resonance lock situations in 
white dwarfs, the three frequencies shown are very close together, and the 
basin of attraction barely extends beyond cj_. 



where /i = MM' j(M +M') is the reduced mass. Since /ia 2 ^> 
/*, in the present context we can neglect any tidal influence 
on Cl even in this extreme-resonance scenario; this also now 
justifies dropping the companion's contribution to E t a e . We 
can then approximate the resonance lock criterion as 



T 
h 



(9) 



This implies that in a resonance lock, the torque increases as 
TixSl 11 ' 3 as the orbit decays. 

Once equation (|9]l is satisfied, the lock can persist for a long 
period of time, since the quantities in equation <|5j affecting 
the magnitude of the torque other than 8uj„ change only very 
gradually in time. (We verify this using numerical orbital evo- 
lution simulations in § [8]) In other words, a resonance lock 
represents a dynamical attractor; Figure|2]provides an illustra- 
tion of the nonlinear dynamics behind a resonance lock. The 
lock could eventually be destroyed if the mode responsible 
began to break near its outer turning point, which would dras- 
tically increase the effective damping rate. This phenomenon 
is explained in §|6] 

If we assume that the system begins completely unsynchro- 
nized, so that a = 2S1, we can determine the orbital period 
where a lock first occurs, which we denote P t \, by substitut- 
ing equations |5) and ([TJ into Q, setting uj n = a and hence 
8uj„ = 0, and then solving for the orbital period. To this end, 
we invoke the following approximate scalings for the eigen- 



TABLE 2 
White dwarf tidal parameters 



ID 


P A (min) 


P trl (min) 


A/ 10- 2 


St 1 100 mill 


He 10 


67 


49 


1.2 


2x 10 9 


He 7 


270 


49 


2.4 


1 x 10' 


He5 


1,400 


90 


4.6 


1 x 10' 


C012 


31 


22 


6.3 


1 x 10 7 


C06 


170 


40 


7.4 


1 x 10 7 



For each of our fiducial white dwarf models from Tab le fT] we list the orbital 
period P A of its first resonance lock (§[4] equation [K)) , the orbital period P tA 
below which traveling wave resonance Tocks can occur (§ |7.2| equation |33) , 
its value of A (§ l5j equ ation |14) , and the tidal quality factor y t | loo min for a 
resonance lock (^ [5.1| equation 1 1 7| evaluated at P OI \, = 100 min. All values 
are for an equal-mass companion. We determined P A and P a \ by directly 
searching over numerically computed eigenmode properties. 

mode linear tidal overlap integral Q and damping rate 7: 

G~ Qo(cr/u*) a and -fRij Q ((T/uj*)~ h , 

where values of the various parameters in these expressions 
are listed in Table|3]for our fiducial white dwarf models. Scal- 
ing parameter values to those for our CO 6 model, we have 



P A ~ 170 min 



2.9 s 



(10) 



where f* = (R 3 /GM) l l 2 is the white dwarf's dynamical time, 
the factor F r \ is 



M'\ fl+M'/M 
~M 



-5/3 



0.010/ \0.i8MR 2 



(id 



70 



2.9 x 10" 1 V 

26^ fn nmn\l/p 



x (2.15 x 10 26 ) (0.00327) l ' p , 
the power p is in general 

P= 473 + L^ <<1 ' (12) 

and p = 0.094 for our CO 6 model (Table B). (The last line of 
equation 1 1 is equal to unity for p = 0.0947) 

For comparison, direct numerical evaluation of eigenmode 
properties with our CO 6 white dwarf model yields P r \ = 170 
min for an equal-mass companion, due to an n = 122 g-mode. 
This is in very good agreement with the analytic approxima- 
tion in equation (JTOj, and is also consistent with our numerical 
results in §[8] We provide values of P r \ for each of our fiducial 
models in Table |2| These results show that resonance locks 
begin at longer orbital periods for cooler white dwarf models, 
due to larger outer convective zones and longer overall dif- 
fusive damping times (smaller j n ; see Table |3]& Figure [7]). 
This increases the maximum possible tidal torque, which is 
proportional to l/7„ (equation ]5V 

5. ENERGETICS 
5.1. Tidal quality factor 
A star's tidal quality factor Q t can be defined as 



Qt = 



£tide 



(13) 
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where the energy c onten t of the tide E^ is approximately 
given by (Appendix B.2 1 



£tide = Ae 



(14) 



A = 2W 2 Y] Q 2 and values of A for various white dwarf mod- 
els are given in Table El Using the relationship bet ween the 
tidal torque and energy transfer rate from equation ( |B8[ ), we 
see that the tidal torque can be expressed in terms of Q t by 



T = 



^tide 



(15) 



Since E^ ~ AF t ide h, where h ~ eP is the height of the equi- 
librium tide and ~ eE* /R is the tidal force, we see that 
our definition of Q t is consistent with 1 /Q t representing an 
effective tidal lag angle; see e.g. |Goldreich & Soter (1966). 
Note that since Qt parameterizes the total tidal energy depo- 
sition rate, which includes mechanical energy transfer asso- 
ciated with increasing the white dwarf spin, the value o f Q t 
alone does not fully determine the tidal heating rate; see § |5.2| 
Using the previous three equations along with the resonance 
lock condition from equation |9j, we have that the value of Q t 
during a resonance lock is 



s.= 



Substituting further yields 



1st 



(16) 



Q t « 9.7 x 10 6 ( -J°*—\ ^ (J*- x ' ' 
1 I 100 min 7 V2.9s, 



M'\ (l+M'/M 



M 



-5/3 



A 

0.074 



(17) 



0.1 8MR 2 



0.010 



Here 7* is the moment of inertia, (3 2 = GM/Rc 2 , and all val- 
ues have been scaled to those appropriate for our CO 6 model 
(Tables [T]& |2) 

Equation ( |17[ ) is a central result of this paper. It is inde- 
pendent of eigenmode properties and is only weakly depen- 
dent on the orbital period, although it depends strongly on 
the mass, radius, and companion mass. Eigenmode properties 
do of course dictate when this value of Q t is applicable, i.e., 
whe n res onance locks are able to occur. We will further show 
in § 7.2 that equation ( fT7] > can hold even when the dynami- 
cal tide is a traveling wave, and the standing wave formalism 
presented thus far is invalid. 

Values of the various quantities entering into equation ( fT7) i 
are provided for a selection of helium and carbon/oxygen 
white dwarf models in Table [T] In particular, since the in- 
spiral time is much longer for low-mass helium white dwarfs 
than for more massive carbon/oxygen white dwarfs, equa- 
tion ( fT7] > predicts that the tidal quality factor Q t should be 
much larger (~ lOOx) for helium white dwarfs, as shown 
in Table [5] meaning tidal effects are more efficient in car- 
bon/oxygen white dwarfs. 

5.2. Tidal heating 

The rate at which heat is dissipated in the white dwarf 
assuming solid-body rotation can be derived using equa- 



tion dB8 



^heat — ^tide ^mech 

d (\ 

dt 



I Q ■ 

2 * spin 



E^Sft 



(18) 



Qt 



where 5fl = Q 

a 



■ n 



spin- 



During a resonance lock we have 



sp m ^iCl = Cl/t gw , so that 

£heat ; 



LCI 5CI 



(19) 



with 5fl then being approximately constant (having negl ected 
rotational modification of white dwarf eigenmodes; see § 9.2 1. 
Defining the asynchronicity period as 5P = 2ir/Sfl, we can 
evaluate this further as 



£he 



L L 

-11/3 



w 

~M 



1.4 x 10 

^orb 

10 min 

h_ 

0.18MR 2 J V0.6M, 



1 +M'/M 



-1/3 



5P 



200 min 



(20) 



M 



8/3 



R 



0.0137? 



o 



again scaling variables to our CO 6 model's properties (Ta- 
ble [T). 

As a simple analytical estimate, consider the example of 
a resonance lock beginning with the white dwarf unsynchro- 
nized at an orbital period Po and continuing until the Roche 
period of PR OC he ~ t * <C Po- The total orbital energy dissipated 
in the white dwarf as heat in this example is 



heat ' 



2nl* 

t 4ir 2 h 
t*Po 

> 7.0E, 



-dt 



(21) 



0A8MR 2 ) 



which could be very large depending on the value of Po- If for 
Po we use our estimate from §|4]of P x \ ~ 170 min appropriate 
for our CO 6 model, we have A£h ea t ~2x 10 47 ergs, a factor 
of ~ 3 larger than the CO 6 model's thermal energy. 

Tidal heating can directly add to a white dwarf's luminos- 
ity and minimally affect its thermal structure if a) the thermal 
time f tn = pc p T/gF at the outer turning point, where wave 

damping is most efficient, is smaller than f gw , and b) £h ea t "C L. 
The outer turning point occurs due to the outer convection 
zone, so ? t h|rcb (radiative-convective boundary) is an appro- 
priate value to use. We find ?th|rcb ^ 10 6 years for all of our 
white dwarf models, as shown in Table [Tj this is -C f gw for 
^orb ^ 10 min, implying that criterion (a)is satisfied. More- 
over, for all models other than our T e ff = 5, 100 K helium and 
J e ff = 5 , 500 K carbon/oxygen white dwarfs, f tn | lc b < 50 years, 
and f tn | rc b < t gw is satisfied even directly prior to mass transfer. 

Criterion (b) above is more restrictive: examining equa- 
tion pO] ) shows that near orbital periods of ~ 10 min, the 
tidal heating rate approaches typical white dwarf luminosi- 
ties. The orbital period where this occurs depends weakly 
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on the various parameters appearing in equation ( |2"0"| ), since 

—11/3 

£heat ^ ^orb • Thus the thermal structure of white dwarfs in 
close binaries may adjust significantly to accommodate the 
additional heat input for f or b < 10 min. We discuss the con- 



sequences of tidal heating further in § 9.2 



5.3. Tidally enhanced orbital decay 

Although the rate f or b at which the orbital period of an 
inspiraling white dwarf binary decays is dominated by the 
gravitational wave term P gw = -P or b/t gv ,, tidal energy dissi- 
pation implies a small deviation from this value (see also Piro 
201 1) ^1 We can compute this difference for a system consist- 
ing oftwo white dwarfs both und ergoing resonance locks by 
using equations |7|, and (jB8]), which yield 



where 



Aide — 



Porb — Pgw Aide ) 



i-s-r gw ' 5w 2 



h+h 
jia 1 



(22) 



(23) 



I\.2 are the moments of inertia of the two white dwarfs, 
[i = M\M2j(M\ +M2) is the reduced mass, a is the semi-major 
axis, and we have again neglected rotational modification of 
white dwarf eigenmodes (i.e. dto n /dfl s ^ = 0, where oj„ is a 
corotating-frame eigenfrequency; see § |9.2[ i. This effect may 
be detectable in future obser vations of close white dwarf bi- 
naries, as discussed in §|9.1| 



6. APPLICABILITY OF STANDING WAVES 

6.1. Wave breaking 

It is important to determine whether the dynamical tide we 
are attempting to study represents a standing wave or a trav- 
eling wave. If it is a standing wave, meaning it is able to 
reflect at its inner and outer turning points without being ab- 
sorbed, then it can achieve large amplitudes due to resonances 
with tidal forcing frequencies (as assumed in §|4}. In the ab- 
sence of nonlinear effects that can occur at large amplitudes, 
a standing wave's damping rat e is w ell approximated by the 
quasiadiabatic value (Appendix |B.3| l, which is small for white 
dwarfs due to their high densities and long thermal times (see 
Figure [7]). On the other hand, if the dynamical tide instead 
behaves as a traveling wave, resulting from absorption prior 
to reflection, then its damping time is approximately a group 
travel time. 

In this section we determine whether the nonlinear process 
of gravity wave breaking causes tidally excited g-modes in 
white dwarfs to be absorbed near the surface and hence to be- 
come traveling waves, as has b een suggested in recent stud- 
ies ( |Fuller & Lai][20TT] [20L2| . Gravity wave breaking has 
been considered extensively in the atmospheric scien ce com- 
munity , since it occurs in Earth's atmosphere; see e.g. |Lindzen| 
( 1981 ). It is also thought to occur in the cores of solar-type 



stars ( [Goodman & Dickson|1998[|Barker & Ogilvie|2010| >. 

Breaking occurs when a wave s amplitude becomes large 
enough to disrupt the stable background stratification. One 
way to derive the condition under which this happens is to 
determine when a wave would produce its own convective in- 
stability, which is equivalent to the perturbed Brunt-Vaisala 

6 Note that the purpose of this section is to determine the influence of the 
tidal energy deposition term Enfe on the rate of orbital decay, even though 
this term is neglected everywhere else in this work, as justified in §|4] 



frequency (squared) becoming comparable to the background 
value — this then makes the total value negative, implying 
convection. The Eulerian perturbation to N 2 is given in linear 
theory by 

^ „A-*-4 + * (24) 



N 2 



where H p is the density scale height and k r is the wavenumber 
in the direction of gravity. Since k r £ r is much larger in mag- 
nitude than the other terms for g-modes, the wave breaking 
condition thus becomes 



(25) 



Other nonlinear processes also come into play when 
\k r £ r \ ~ 1. Indeed, this criterion is equivalent to Ri ~ 1/4, 
where Ri is the Richardson number due to the wave's shear, 
which implies the wave is Kelvin-Helmholtz unstable. Equa- 
tion ( p5j ) is also similar to the condition under which surface 
ocean waves break: when the vertical displacement becomes 
comparable to the wavelength. 

To determine whether g-modes break, we evaluated the lin- 
ear, quadrupolar tidal fluid response assuming global adia- 
batic normal modes and an equal-mass companion; see Ap- 
pendix [B] Under these assumptions, we find that for both our 
helium and carbon/oxygen white dwarf models, the dynami- 
cal tide breaks for close resonances at orbital periods as large 
as ~ 1 hr, as shown in Figure|3] The off-resonance dynamical 
tide begins to break more generically at P OI b < 10-20 minQ 
Furthermore, for all of the white dwarf models we have con- 
sidered (Table[TJi, we find that at sufficiently long orbital peri- 
ods, the dynamical tide doesn't break even for a perfect reso- 
nance, and thus that standing wave resonance locks should be 
able to occur. As such, we expect wave breaking not to oper- 
ate during a significant portion of the inspiral epoch, in which 
case the analysis presented in §|4]may be valid. We address 
this in more detail in §[8] 

6.2. Differential rotation and critical layers 

The possibility of differential rotation represents a signif- 
icant challenge to the standing wave assumption we utilized 
in § [4] Indeed, tidal angular momentum is preferentially de- 
posited in the outer layers of a white dwarf, since that is where 
damping times are shortest and waves are able to communi- 
cate their energy and angular mo mentum content t o the back- 
ground stellar profile ( |Goldreich & Nicholson| 1989 b). Thus 
tides do not naturally induce solid-body rotation, and instead 
tend to first synchronize layer s near the outer part of the grav - 
ity wave propagation cavity ( |Goldreich & Nicholson| 1 989a), 
absent the influence of efficient internal angular momentum 
transport. 

The presence of a synchronized or "critical" layer at the 
edge of a mode propagation cavity implies that the mode's 
coro taring frequency tends to zero at that location, which in 
turn means its radial wavenumber becomes very large due 
to the asymptotic g-mode dispersion relation cj ~ N(kh/k r ), 
where k^ and k r are respectively the perpendicular and radial 



7 This is in conflict with the claims made in |Fuller & Lai|(2012} , since 
that work used k r \£\ ~ 1 to assess wave breaking, instead of :quation {25}. 
The total displacement |£| = + £jj)'/ 2 includes horizontal motion, which 
is perpendicular to the stratification and thus does not contribute to breaking. 
As a g-mode's horizontal motion is much greater than its vertical motion, they 
overestimated the degree of breaking by a factor of ~ £h/£r ~ 3> 1, 

where a = 2(Q- Q s „- m ) is the / = m = 2 tidal driving frequency. 
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FIG. 3. — The maximum value of \k r £ r \ attained throughout the propagation cavity of our C06 (black line) and He5 (red line) white dwarf models (Table^, 
assuming adiabatic standing waves, f2 S pin = £1/2, and an equal-mass companion. Where |<r r £r| m ax > 1, wave breaking occurs (§ |6.1| , and the effective wave 
damping time becomes roughly the group travel time across the white dwarf. This occurs very near resonances for P or (, < 0.5 - 1 hr. 



wavenumbers. As a result, the mode's local damping time 
becomes very short, and it is absorbed rather than reflected, 
eliminating the possibility of achieving resonant amplitudes 
(although traveling waves can also effect resonance locks at 
short orbital periods; §[7J. 

In Appendix [A] we analyze angular momentum redistribu- 
tion by fossil magnetic fields, possibly generated by a progen- 
itor star's convective core (during hydrogen or helium fusion) 
and amplified by flux freezing as the core contracts. We calcu- 
late that a field strength of only ~ 200 G is required to main- 
tain solid-body rotation during a resonance lock for an orbital 
period of ~ 100 min in our CO 6 mod el, and onl y ~ 20 G 
in our He 7 model (Table [TJ. |Liebert et al.| ( |200"3) conclude 
that at least ~ 10% of white dwarfs have fields > 10 6 G, and 
speculate that this fraction could be substantially higher; field 
strengths in white dwarf interiors may be even more signif- 
icant. With a field of 10 6 G, our calculations indicate that 
critical layers should not occur until orbital periods of less 
than 1 min, or even less if the field can wind up significantly 
without becoming unstable. 

6.3. Validity of the secular approximation 

The Lorentzian mode amplitude solutions invoked in Ap- 
pendixjBjto produce the standing wave torque in equation |4j 
are strictly valid only when a mode's amplitude changes 
slowly relative to its damping time. Further examining equa- 
tion |5j, we see that near a perfect resonance the amplitude 
changes by a factor of ~ 2 as the detuning frequency changes 
by of order the damping rate j n . Thus the Lorentzian solution 
is applicable near a perfect resonance only when 



t" 1 < t — 

In ~ 'gw o ' 



In 



which evaluates to 



fU > 90 min 



M \ 



5/11 



O.6M0 



(26) 



(27) 



where 



M 
~M 



3/11 



1 +M'/M 



-l/ii 



In 

80 yr 



6/11 



Equation ( |27[ i is scaled to values for our CO 6 model (Fig- 
ure pj; the restriction instead evaluates to P 01 -b > 50 min for 
our He 7 model, using a damping time of j' 1 ~ 60 yr ap- 
propriate for the initial resonance lock. Below these periods, 
the Lorentzian solution becomes invalid and the exact out- 
come is unclear, although our preliminary numerical integra- 
tions of fully coupled mode amplitude and orbital evolution 



equations indicate that resonance locks can still occur even 
beyond the validity of the Lorentzian solution. (We address 
a similar conce rn relating to angular momentum transport in 
Appendix A. 2 ) Nonetheless, we find that the initial standing 
wave resonance lock occurs at orbital periods larger than the 
critical value from equation ( |27"| ) in our CO 6 and He 7 mod- 
els (Table [5J, meaning resonance locks should proceed as ex- 
pected. 

7. TRAVELING WAVES 

7.1. Excitation and interference 

In this section, we will describe two different mechanisms 
of tidal gravity wave excitation considered in the literature. 
We will then compare both sets of theoretical predictions to 
our numerical results to assess which mechanism predomi- 
nantly operates in our fiducial white dwarf models. 



Zahn| ( |1975| l showed that when a gravity wave is well de- 
scribed by its WKB solution, a conserved wave energy flux 
results. Thus gravity waves must be excited where the WKB 
approximation is invalid: where the background stellar model 
— particularly the Brunt- VaisaTa frequency N — changes 
rapidly relative to a wavelength. 

One natural candidate for wave excitation, then, is at a 
radiative-convective boundary (RCB), where N 2 abruptly be- 
comes negative. White dwarfs possess convective envelopes 
near their surfaces (Figure [TJ, so thi s mechanism is plausible. 
The res ulting theoretical prediction ( Za hn|197"5| |Goodman & 
Dickson 1998) is that the traveling wave tidal torque should 



scale as r oc c 8 / 3 , where a = 2(f2-J7 sp i n ) is the m = 2 tidal 
driving frequency. Using our c alcu lation of the traveling wave 
torque in equation (|32|i from §|7.2| we see that th is in turn im- 



plies that the linear overlap integral (Appendix B.4i should 

scale as Q„ o c u l „ X ^ ', given our normalization convention in 
equation ( |B5| l. 

More recently, [Fuller & Lai ( 2012} showed that excitation 
can also proceed near the spike in the Brunt- VaisaTa frequency 
that occurs at the transition between carbon/oxygen and he- 
lium in a carbon/oxygen white dwarf (see Figure [TJ. Their 
corresponding prediction for the torque scaling is t oc a 5 , im- 
plying Q n oc w 3 . Thus this mechanism predicts a steeper over- 
lap scaling with frequency than for excitation at the RCB. 

An additional feature of excitation at a composition bound- 
ary is that waves originate from a location inside the propaga- 
tion cavity, meaning both an ingoing and outgoing wave are 
created. Since the ingoing wave reflects at the inner turning 
point, interference occurs between the reflected ingoing wave 
and the purely outgoing wave. Constructive interference im- 
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FIG. 4. — Plots of the linear tidal overlap integral Q„, which characterizes the spatial coupling strength between the tidal potential and a given mode (Ap- 
pendix |B~4) , as a function of the eigenmode frequency ui„ and radial order n, for the first 500 g-modes in four of our fiducial white dwarf models (Table]!}. Panels 

1-3 are helium white dwarfs ordered by increasing temperature; panel 4 is a carbon/oxygen white dwarf. A smooth power law scaling of Q„ oc oj„ im plies 
that gravity wave excitation by the tidal potential occurs at the interface between a white dwarf's outer convection zone and its inner radiative core (§ |7.1) ; this 
can be seen in the cooler h elium models from panels 1 & 2. Hotter white dwarfs have smaller convective regions, and wave excitation instead may occur at 
composition gradient zones ( Fuller & Lai 2012 1; this mechanism predicts steeper, more jagged profiles of Q„ with ui n , as in panels 3 & 4. 



plies a large overlap integral Q n , whereas destructive interfer- 
ence makes the overlap small, thus this mechanism predicts a 
jagged overlap profile with respect to the wave frequency uj„. 

With these theoretical predictions in hand, the essential 
question to answer is which excitation mechanism — RCB 
or composition gradient — is most efficient in a given white 
dwarf modelj^] The answer hinges on the properties of the 
convective envelope. Table [T] and Figure [T] show that this en- 
velope is very small in hot white dwarfs, and exists at very low 
densities, but that its extent increases rapidly as a white dwarf 
cools. Thus it seems possible that excitation at the RCB may 
occur for cooler white dwarfs, whereas hotter white dwarfs 
must rely on the composition gradient mechanism. 

One method we can utilize to distinguish between the 
two mechanisms is simply to observe the power law scal- 
ing Q n oc oj° of numerically computed linear overlap inte- 
grals for our various white dwarf models, given in Table [3] 
Consistent with our expectations, cooler helium white dwarfs 
have a power law index a ~ 1 .83 « 11/6, implying excitation 



TABLE 3 

White dwarf I = 2 eigenmode properties 



ID 


Qo 


a 


7o/u>* 




b 


cco/u* 


c 


HelOl 


9.6 x lfr 6 


2.61 


1.5 x 10- 


] i 


6.16 


0.0891 


2.00 


He7 


3.6 x IO" 6 


1.83 


2.1 x 10" 


12 


2.00 


0.158 


2.00 


He5 


7.8 x 10~ 4 


1.90 


7.7 x 10" 


15 


1.99 


0.298 


2.00 


C012t 


7.2 x 10 1 


4.40 


1.2 x 10" 


14 


6.41 


0.403 


2.00 


C06 


2.7 X 10 1 


3.69 


2.9 x 1(T 


14 


1.88 


0.743 


2.00 



8 This question was not addressed in |Fuller & Lai| j2012| since they 
adopted an absorbing boundary condition near the outer turning point, and 
thus did not include the convection zone in their calculations. 



Asymptotic fits to numerically computed eigenmode properties for the white 
dwarf models from Table [T] The linear overlap integral Q„ (Appendix B.4 
& Figurepl is fit as Q„ =Qo(u) n /u)*) a ; the damping rate j a (Appendix B3 
& Figurel7l is fit as 7,, = 7o(i^ , 1 /a; t )~ fc ; and the inverse group travel time 
a„ = 27r/fg rol ip „ (Appendix B.3| is fit as a„ = ao(^ji/w H ,) £ '. tNote that rapid 
thermal diffusion near the outer turning point causes g-modes of radial order 
n > 50 to become traveling waves in our C012 and He 1 mod els, meanin 
our fits for 7,, are not relevant in this regime; see Appendix|B.3|and Figuref 



at the RCB, while hotter helium white dwarfs and our car- 
bon/oxygen models have larger values of a. Furthermore, Fig- 
ure]?] shows that models with steeper overlap power laws also 
show jagged variation of Q n with frequency, thus demonstrat- 
ing the interference predicted by composition gradient excita- 
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tion. 

Note that if a gravity wave in a cool heliu m wh ite dwarf be- 
gins to break near its outer turning point (§ 6.1 1, this implies 
that the tidal excitation and wave breaking regions would be 
almost directly adjacent. It might then be possible for break- 
ing to inhibit excitation, meaning that the composition gradi- 
ent mechanism would again dominate. A more sophisticated 
hydrodynamical calculation is required to address this con- 
cern. 

7.2. Traveling wave resonance locks 

The resonance lock scenario we described in §|4]relied on 
resonances between standing white dwarf eigenmodes and the 
tidal driving frequency. However, resonance locks are in fact 
a more general phenomenon that does not explicitly require 
standing waves P] In the context of white dwarf binary inspiral, 
the two essential requirements on the tidal torque function r 
in order for a resonance lock to occur are: 

a) The torque profile must be a jagged function of the 
l = m = 2 tidal driving frequency a = 2(fi-r2 sp i n ) (equa- 
tion[3"T|below). 

b) The magnitude of the tidal torque must be large enough 
that it can satisfy equation |[9): t = hVl/t^. 

When these conditions are satisfied and a resonance lock oc- 
curs, the tidal quality factor Q t and heating rate are given by 
equations ( fPTj ) and pO) , respectively. 

We first address criterion (a). Dropping the tidal energy 
deposition term from equation |7]), as justified in § [4] yields 
the simplified orbital evolution equation 



1 da 
m dtt 



(28) 



where the gravitational wave decay time f g „(ST) is defined in 
equation ([!}. 

Let us assume that the tidal torque satisfies criterion (b) at 
an orbital frequency f^o and a tidal driving frequency <jq, so 
that da/dil = and equation ( |28] i reduces to 

4f>o = fgw(^oMr!o,CTo). (29) 
As long as r increases with a, equation d29l) represents a sta 



ble fixed point of the evolution equations; see Figure [2] Next, 
since the orbital frequency steadily increases due to the emis- 
sion of gravitational waves, we examine what happens to this 
fixed point when changes by a small amount +Af2. In or- 
der to preserve equation d29j, the tidal driving frequency must 
commensurately change by an amount Aa given by 



Act _ 1 (An\ /dlogr 
~a~ ~~3 \dloga 



(30) 



which can be derived by differentiating equation ( |29] l and sub- 
stituting equation Q. 

Equation ( |30| ) allows us to appropriately quantify the 
"jagged" variation of the torque function required by criterion 
(a): if 

dlogr 



> 1, 



(31) 



d loger 

9 "Synchronization lock" may be a better term in the traveling wave limit; 
however, we continue using "resonance lock" due to the many similarities 
between standing wave and traveling wave results. 



then the fixed point can be maintained by only a minimal 
change in the forcing frequency for a given increase in the 
orbital frequency, thus constituting a resonance lock. Any 
general power law trend of r with a will fail to satisfy this 
condition — additional sharp features are required^] 

Torque profiles consistent with equation ( [3T] > can be pro- 
vided in several ways. For standing waves, the comb of 
Lorentzians produced by resonances with eigenmodes (see 
equation |4] and Figure [3]) easily satisfies equation ( [31) , since 
white dwarf eigenmodes are weakly damped, meaning on- 
and off-resonance torque values differ by many orders of 
magnitude. For t raveling waves, if th e composition gradi- 
ent mechanism of [Fuller & Lai| ( 2012) discussed in § |7.1| is 
the dominant source of wave excitation, it naturally provides 
sharp features in the torque function due to wave interference. 
This can also be observed in Figure [5] where the traveling 
wave torque changes by a factor of ^5 as a = 25Q changes 
by only ~ 10%, implying |c/logT/t/logcr| ~ 50. 

Lastly, wave breaking can also provide rapid variation in 
the torque profile due to a sudden transition between stand- 
ing and traveling wave torques that occurs near resonances 
at short orbital periods. Specifically, as the tidal driving fre- 
quency a sweeps towards a resonance due to orbital decay 
by gravitational waves, a tidally excited g-mode's ampl itud e 
can become large enough to induce wave breaking (§ 6.1 1, 
which causes the effective damping rate and hence the result- 
ing torque to increase enormously (see the blue curve in Fig- 
ure |5j. The precise shape of this transition requires hydrody- 
namical simulations to ascertain; fortunately, we find that es- 
sentially any transition between a nonresonant standing wave 
torque in between resonances and a traveling wave torque near 
resonance will satisfy equation pT) for white dwarfs, due to 
the large disparity between typical damping times associated 
with standing waves and the group travel time, which approx- 
imates the damping time for a traveling wave (Table [3] & Fig- 
ure|7]i. We discuss this further in §[8] 

Next, we address criterion (b) for a resonance lock stated 
at the beginning of this section by estimating the mag nitude 
of the traveling wave torque. Following Goodman & Dickson 
( 1998 ), we employ the standing wave result for the torque due 
to a specific mode from equation |5]l with the tidal driving 
frequency a set to a particular eigenfrequency uj n , j n replaced 
by a„ = 2 7r/f gr oups (where f gl - oup is the group travel time; see 
Appendix B.3 1, and 6lj„ set to zero. This yields 



T tt av(o- = w„ , fi) ~ 4E* e Q n uj n / a n 



(32) 



where we have approximated W 2 ~ 1 . Then, in order to evalu- 
ate an effective traveling wave torque for arbitrary a, we sim- 
ply interpolate over values computed using equation ( j32[ ). 

To estimate the first orbital period P v \ at which traveling 
wave resonance locks can occur, we follow the same proce- 
dure as in § |4] and again invoke approximate scalings for the 
eigenmode linear tidal overlap integral Q (Appendix B.4) an d 
the effective traveling wave damping rate a (Appendix B.3 1: 



Q ~ Go(o-/w*) fl and a w a (a/u*T, 



10 [Fuller & Lai|f2012[ also noticed that a constant occurred in their 
simulations, although they attributed this to the overall power law trend of 
their torque function with a. Indeed, their results possess sharp interference- 
generated features that provide a much larger contribution to |rflogT/rflogcr| 
than the trend, meaning a resonance lock was likely responsible for maintain- 
ing <r pb constant. 
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where c = 2; see Tableland Figured The resulting formula, 
scaled to values for our CO 6 model (Table [TJ, is 



' 43 min 



(—) 
\2.9s) 



(33) 



where f* = (R 3 /GM) X I 2 is the white dwarf's dynamical time, 
the factor F tr \ is 



1 +M'/M 



-5/3 



0.010/ \0.18MR 2 

2 / , 
ceo 



(34) 



Go 

27 J V0.74w* 
(8.41 x 10 12 ) (0.0119) 17 



the power q is in general 



q = 



1 



-1/3 + 2a 



< 1, 



(35) 



and q = 0.15 for our CO 6 model (Table BJ. (The last line of 
equation [34] is equal to unity for q = 0.15.) Equation ( [33) 
assumes the white dwarf begins c ompletely unsynchr onized; 
it is equivalent to equation (79) of [Fuller & La i (2012). 

Direct numerical evaluation of eigenmode properties with 
our CO 6 white dwarf model yields P a \ = 40 min for an equal- 
mass companion, due to an n = 27 g-mode, which agrees well 
with equation ( |33) . Values of P tl -i for each of our fiducial mod- 
els are provided in Table [2] Note, however, that in deriving 
these results for f tr i we have assumed that the white dwarf 
spin is much smaller than the orbital frequency; if significant 
synchronization has already occurred, the true value of P a \ 
will deviate from our prediction by an order-unity factor. 

8. NUMERICAL SIMULATIONS 

To address the tidal evolution of an inspiraling white dwarf 
binary undergoing resonance locks, we aim to combine the 
standing and traveling wave results from §§ [4] & [71 numeri- 
cally. To this end, we evaluate the complete standing wave 
tidal torque from equation |4) and solve for the spin and 
orbital evolution using equation (JTJi. To account for wave 
breaking, we check that all ei genm odes satisfy |fc r 6-|max < 1 
throughout the white dwarf (§ |6.1| i; when an eigenmode ex- 
ceeds unit shear, we inst ead set its damping rate to a„ = 
27r/f g roup,n (Appendix B.3 i, which approximates the traveling 
wave regime (Goodman & Dickson 1998). We smoothly tran- 
sition between the standing and traveling wave regimes using 
the interpolation formula 



T = 



: 

1 + ( |&r£r|max) 



(36) 



where r stan d is the standing wave torque from equation (|4), 
is the traveling wave torque produced by interpolating 
over equation ( [32) , and | | max is the maximum value of 
the wave shear over all rele vant eigenmodes and across the 
entire propagation cavity (§ |6.1[ ), evaluated assuming stand- 
ing waves. We arbitrarily adopt z = 25 to induce a sharp 
transition that occurs only when |£ r £,.| max is very close to 1; 
our results are insensitive to the value of z so long as it is 



lO" 12 
£ lO" 14 
§ 10" 16 

r io- 18 

| IO" 20 
10" 22 

50.0 

10.0 
1 5.0 




1.0 
0.5 

0.1 



stahding 



.traveling, 
■transition ■ 




0.23 0.24 0.25 0.26 
6fl ( rad / min) 



0.27 



FIG. 5. — Top panel: Example plot of the standing wave torque r stan j (equa- 
tion HI red line), the traveling wave torque T tl - av (interpolation over equa- 
tion R5] evaluated at eigenmode frequencies; d ashe d green line), and our in- 
terpolation between the two regimes (equation |36| thick blue line), as func- 
tions of <5f! = CI — f2 sp j n at fixed P or |, = 30 min tor our C06 model (Table^ 
and an equal-mass companion. The standing wave torque on average is many 
orders of magnitude smaller than the traveling wave torque; however, near 
resonances it becomes many orders of magnitude larger. Wave breaking acts 
to "cap" the Lorentzian peaks of the standing wave torque in the interpolation 
function. Bottom panel: Plot of the wave breaking criterion \k r ^ r \ maximized 
over all eigenmodes and the entire propagation cavity (blue line), using the 
same parameters and model as the top panel. When | | max < \ t the dy- 
namical tide represents a traveling wave, and the tor que r — > T stail d; when 
|fcr£r|max > 1> wave breaking occurs, and r — > r trav (§ |6.1| . 

> |ln(r trav /r st and)|- Figure^shows a comparison of r stand , r trav , 
and the transition function in equation ( [36) . 

Figure [6] shows the results of two of our simulations. The 
left column used our O.2M , T e fs = 7, 000 K He 7 model, while 
the right column used our O.6M , T e ff = 5,500 K CO 6 model 
(Table[T). We did not account for white dwarf cooling or tidal 
heating, and instead used fixed white dwarf models through- 
out both simulations. We initialized our simulations with 
^spin ~ 0, and the orbital period set so that the time until mass 
transfer f merge = 3f gw /8 (equation [T]i was equal to 10 billion 
years. 

Both simulations follow the archetypal scenario laid out in 
§[3] transitioning amongst the four regimes (SI), (S2), (Tl), 
and (T2). Both begin in (SI), where the dynamical tide is a 
standing wave even near resonances, but the tidal torque is 
too weak to create a resonance lock. As the orbit shrinks due 
to gravitational wave radiation, the tidal force waxes and the 
first resonance lock eventually begins in both simulations at 
the appropriate value of P r \ estimated in § [4] and provided in 
Table [2j this is regime (S2). At this point tidal heating and 
synchronization suddenly become much more efficient (§[5), 
and the difference between orbital and spin frequencies re- 
mains constant. Analytic formulas for the tidal quality factor 
and heating rate appropriate for this situation (as well as T2 
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He 7 white dwarf model: T eS = l, 000 K, O.2M CO 6 white dwarf model: T eft = 5 , 500 K, 0.6M, 
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FIG. 6. — Results of numerical simulations of the secular evolution of white dwarf binaries, the details of which are described in §[8] The left column shows 
results using our 0.2Mq, T t g = 7,000 K He7 model, while the right column used our 0.6Mq, r e ff = 5,500 K C06 model (TablefT). Top row: The orbital 
(dashed black line) and spin (thick blue line) frequencies, as well as their difference 8Q = f2 — f2 s pi„ (red line); the latter sets the tidal forcing frequency a — 28Q. 
Resonance locks correspond to regions where 8Q is constant. Our assumption of slow rotation breaks down when 5f2 < C sp j n due to nonlinear rotational 
modification of stellar eigenmodes. Second row: Number of radial nodes n of dominant eigenmo de/w ave. Third row: Maximum value of \k r £ r \ across entire 
white dwarf, evaluated assuming standing waves, which assesses whether wave breaking occurs (§ |6.1| . During the initial resonance lock, |^ r £r|max starts < 1, 
but gradually rises until it becomes ~ 1 and breaking begins. Fourth row: Rate at which orbital energv is dissipated as heat in the white dwar f, in units of Lq. 
Bottom row: Tidal quality factor Q t (blue line) and time until mass transfer f mc rge = 3fg w /8 (equationTTl dashed magenta line). See equations j!7| and < |20| for 
analytic estimates of the tidal quality factor Q t and heating rate, respectively. 
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discussed below) are given in equations ( fT7| i and ( |20| ), respec- 
tively, and exactly reproduce their numerically derived values 
appearing in the bottom two rows of Figure [6] 

Both simulations begin the standing wave resonance lock 
regime (S2) with a value of the wave breaking criterion 
|&r£r|max < 1; however, as the orbit shrinks further, progres- 
sively larger and larger wave amplitudes become necessary to 
support a resonance lock, eve ntual ly leading to wave breaking 
near the outer turning point (§ 6.1 1. At the onset of (S2) in the 
C06 simulation, the inequality |K r £ r | max < 1 is only weakly 
satisfied, meaning that the standing wave lock regime (S2) is 
short lived. In the He 7 simulation, however, (S2) begins with 
l^rCrlmax "C 1, so that the initial resonance lock persists from 
P r \ = 270 min to a period of P m 50 min, corresponding to an 
interval of time of about 5 billion years. 

Once |fc r £r|max becomes ~ 1, both simulations enter regime 
(Tl), where near would-be resonances the dynamical tide be- 
comes a traveling wave too weak to create a resonance lock. 
The otherwise steeply peaked standing wave torque is thus 
capped in this regime; see Figure [5] Regime (Tl) results in a 
weak tidal synchronization and heating scenario, very similar 
to (SI). 

Eventually, at an orbital period ~ P tl i (§ |7.2| Table|2]i, both 
simulations enter regime (T2), where even the traveling wave 
torque can create a resonance lock. Tides again become ef- 
ficient, with synchronization and heating scenarios quantita- 
tively consistent with the analytic results in § [5] (just as in 
S2). In the He 7 simulation, (T2) begins at an orbital period of 
w 27 min, which differs from its value of P tt i = 49 min listed 
in Table [2] since that value is only strictly applicable when 
^spin = 0, whereas significant synchronization has already oc- 
curred. The value of P a \ in the CO 6 simulation is a better 
estimate of the onset of (T2) due to the brief duration of (S2) 
in that case. 

The maximum wave shear |& r £ r |max shown in Figure [6] 
(which is evaluated assuming standing waves) remains very 
close to unity throughout much of regime (T2). A reasonable 
question, then, is whether this is an artifact of the interpolation 
function we used to transition between standing and traveling 
waves torques (equation [36}. 

On the contrary, we believe there is a physical reason why 
|£,£, | max should saturate at ~ 1, and that it is a natural con- 
sequence of t he tr aveling wave resonance lock scenario we 
proposed in § 7.2 Specifically, at this point in the system's 
evolution, if the dynamical tide attempts to set up a standing 
wave, the orbital frequency will evolve, increasing the tidal 
driving frequency a = 2(£1 - £l sp i n ) towards a resonance and 
inducing wave breaking. However, fully transitioning to the 
traveling wave regime then creates a much larger torque (due 
to the much larger effective damping rate), causing the spin 
frequency to increase rapidly and sending a away from res- 
onance, ending wave breaking and reinstituting the standing 
wave regime. The end result is that |fc r £r|max should average 
to be ~ 1. 

This line of reasoning suggests that the true phenomenon 
may be episodic in nature. Alternatively, a weak-breaking 
regime may be possible, allowing the system to smoothly skirt 
the boundary between linear and nonlinear fluid dynamics. 
Full hydrodynamical simulations may be necessary to under- 
stand this in more detail. 

9. DISCUSSION 
9.1. Observational constraints 



The theoretical results we have developed can be 
compared to the recently discovered system SDSS 
J065133.33+284423.3 (henceforth J0651), which con- 
sists of a r e ff = 16,500 K, O.26M helium white dwarf in a 
13-minute eclipsing binary with a = 8, 700 K, 0.50M Q 
carbon/oxygen white dwarf ( [Brown et al. 201 1 1. Orbital 
decay in this system consistent with the general relativistic 



prediction was discovered by Hermes et al. ( 2012| l. 

|Piro| ( |20TT| l studied tidal interactions in J0651, and pro- 
duced lower limits on values of the tidal quality factor Q t by 
assuming that the observed luminosity of each white dwarf is 
generated entirely by tidal heating, and that both white dwarfs 
are nonrotating. However, the definition of Q t used in that 



work differs with ours (equation 13 1, which hinders straight- 
forward comparison^] 

Instead, we compare the observed luminosities with our ex- 
pression for £heat f rom equation ( |20] i, which is applicable dur- 
ing a resonance lock. All parameters for J065 1 entering into 
( f2"0"] i were determined observationally except the moments of 
inertia and the asynchronicity periods SP = 2-k/SQ,, where 
<5S1 = f2-S! sp j n . We can thus use appropriate values of from 
Tableland impose the inequality L > E^ e . dt , since cooling can 
also contribute to each luminosity, in order to constrain SP for 
each white dwarf. 

This calculation yields 8P > 7 min for the helium white 
dwarf and SP > 400 min for the carbon/oxygen white dwarf, 
each with uncertainties of ~ 20%. Since the orbital period of 
J0651 is such that resonance locks should currently exist in 
both white dwarfs — P 01 -b is le ss th an both P r \ from § |4] for 
standing waves and P tI \ from § 7.2 for traveling waves (Ta- 
ble 12} — and since our simulations developed wave breaking 
longbefore P 01 -b = 13 min (Figure |6j, our a priori expectation 
is that SP should be ~ P a \ ~ 50 min. It is encouraging that the 
inferred constraints on SP for both white dwarfs are within an 
order of magnitude of this prediction. 

We can nonetheless comment on what the deviations from 
our predictions may imply. First, the fact that SP > 7 min < 
^orb = 13 min for the helium white dwarf implies that much 
of its luminosity must be generated by standard white dwarf 
cooling or residual nuclear burning rather than tidal heating. 
If this is correct, it would imp ly an age for the helium white 
dwarf of only ~ 40 Myr (Panei et al. 2007); dividing this age 
by a cooling time of ~ 1 Gyr (Table 11} yields a very rough 
probability for finding such a system of ~ 4%. This scenario 
does not seem unlikely, however, since selection bias favors 
younger white dwarfs. 

On the other hand, the lower limit of ~ 400 min placed 
observationally on SP for the carbon/oxygen white dwarf is 
much larger than our predictions for both P r i and P u -i. Further- 
more, tidal heat is deposited very close to the photosphere in 
hot carbon/oxygen white dwarf models (Table [TJ, so we ex- 
pect tidal heating to contribute directly to the l uminosity of 
the carbon/oxygen white dwarf in J0651 (§ 5.2 1. The obser- 



vational constraint on SP thus means that the carbon/oxygen 
white dwarf in J065 1 is less luminous than our analysis pre- 
dicts, by a factor of ~ Pz/Pu\ ~ 10. 

Although this is formally inconsistent with our results, ex- 
amining the first row of Figure [6] shows that both of our nu- 
merical simulations have SQ <C "spin near ^orb =10 min. This 

" The relationship between our value of the tidal quality factor, Q t , and 
that used inlPirol 4201 1|, Q[, is Q[ = QtcrM/XQ/i, where fj, is the reduced 
mass. Our value y t is consistent with being the recipr ocal of an effective 
tidal lag angle, which is the conventional definition; see § |5.1| 
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means that the influence of rotation on eigenmode properties 
is like ly to be very important at such short orbital periods 
(§ 9.2 1, which is not included in our analysis. Damping and 
excitation of white dwarf eigenmodes by nonlinear processes 
are also likely to be important considerations. 

Lastly, assuming resonance locks are occurring in both 
white dwarfs, we predict that t he r ate of orbital decay should 
be enhanced due to tides by (§ 5.3 i 



tide 



g" ' J0651 



2%. 



where f gw = -P rbAgw Although this estimate fails to include 
the effect of rotation on eigenmode frequencies, which we al- 
ready argued may be important in J065 1 , it should nonethe- 
less be robust at the order-of-magnitude level. This ~ 2% 
deviation between the system's period derivative and the gen- 
eral relativistic prediction not accounting for tides may be de- 
tectable given further sustained observations (Pir o|20lT) . 

9.2. Rotation and white dwarf evolution 

In our analysis we have neglected the influence of rota- 
tion on the stellar eigenmodes beyond the simple geometri- 
cal Doppler shift of the forcing frequency into the corotating 
frame. To linear order in the rotation frequency, the correction 
to the stellar eigenfrequencies makes very little difference to 
the results we have derived — it just means there should be 
factors of (1 — C„) ~ 5/6 appearing in various formulas in §|4j 
which we neglected for simplicity. 

However, when 5£l = CI— O sp i n < S! sp j n , nonlinear rotational 
effects become important. Figure |6] shows that this inequality 
is satisfied below f or b ~ 25 min in our CO 6 simulation, and 
takes hold soon after the first resonance lock in our He 7 sim- 
ulation, at only P or b ~ 150 min. Below these orbital periods, 
fully accounting for the Coriolis force in the stellar oscillation 
equations becomes necessary. 

For example, excitation of rotationally supported modes — 
Rossby waves and inertial waves — could prove very effi- 
cient. Such modes have corotating-frame eigenfrequencies 
that are strongly dependent on the rotation frequency, so a res- 
onance lock wou ld follow the more general trajectory fWitte| 
& Savonije|1999l > 



= SUJ: 



mcKV I spm 



L spin 



ci 



(37) 



Since our analysis in § [5] relied on resonance locks produc- 
ing CI m Cl spm , which no longer holds when dLu n /dCt sp i n it 
is unclear whether nonlinear rotational effects could substan- 
tially alter our results for e.g. the tidal quality factor (equa- 
tion 17 1 and tidal heating rate (equation 20 1. 

WnTte dwarf cooling and tidal heating could also poten- 
tially modify the synchronization trajectory that results during 
a resonance lock. For example, since g-mode frequencies ap- 
proximately satisfy oj„i ~ (N)l/n, and since the Brunt- Vaisala 
frequency scales with temperature as N oc T x l 2 in a degener- 
ate environment (§|2|, progressive changes in a white dwarf s 
thermal structure due to either heating or cooling would in- 
troduce an additional doj„/dt term on the right-hand side of 
equation ( f37| ). 

9.3. Crystallization 



Whether a plasma begins to crystallize due to ion-ion elec- 
tromagnetic interactions is determined by the Coulomb in- 
teraction parameter T, which is defined as the ratio of the 
Coulomb to thermal energy, 



2, ,2 



Z z e 
dJcT' 



(38) 



where Ze is the mean ion charge and d\ is the ion separa- 
tion, defined by 1 = rii(4Tr/3)df . When V > 1, the plasma 
under consideration behaves as a liquid; when V > r crys , the 
plasma crystallizes. This critical value is r crvs ~ 175 in single- 
component plasmas (Dewitt et al. 2001 ). However, more re- 
cent observational studies of c arbon/oxygen white dwarf pop- 
ulations (Wing eTet al.||20 09) as well as detailed theoretical 
simulations (Horowitz et al. 2007) indicate that a larger value 
of r crys ~ 220 is applicable for two-component plasmas, as in 
the cores of carbon/oxygen white dwarfs. 

As shown in Table [T] the central value of T does not ex- 
ceed the appropriate value of r crys for any of our helium white 
dwarfs. However, for C06, our 0.6M Q , r eff = 5,500 K car- 
bon/oxygen white dwarf, we have r core = 260 > r crys , and 
further T > r crys for 19% of the model by mass (taking 
r c rys = 220). This is indicated in the bottom panel of Figure[TJ 
as a shaded region. 

The excitation of dynamical tides in white dwarfs possess- 
ing crystalline cores is an interesting problem that deserves 
further study. We will only speculate here on the possible 
physical picture. Our preliminary calculations of wave prop- 
agation inside the crystalline core, using expressions for the 



shear modulus of a Coulomb crystal from|Hansen & van Horn 



(1979), indicate that the shear wave Lamb frequency is sev- 
eral orders of magnitude too large to allow gravity waves to 
propagate as shear waves in the core. 

Thus it seems possible that dynamical tides could be ef- 
ficiently excited at the edge of the core, as in excitation 
at the edge of a convective core in early-type stars ( |Zahn| 
|1975| l. In this scenario the deviation of the tidal response 
inside the crystal from the potential-filling equilibrium tide 
solution excites outward-propagating g-modes. As a conse- 
quence, tidal gravity waves may be much more efficiently ex- 
cited in crystalline core carbon/oxygen white dwarfs, since 
the Brunt-Vaisala frequency gradient near the core would be 
much steep er th an in any composition gradient zone, as dis- 
cussed in § |7.1| 

9.4. Nonlinear damping 

Since eigenmodes responsible for resonance locks in the 
standing wave regime attain large amplitudes, it is natural to 
worry that they might be unstable to global nonlinear damp- 
ing by the parametric instability, even i f they don't experienc e 
wave breaking (e.g. Arra s et al.|2003 Weinberg et al.|20T2] |. 
To address this, we performed a rough estimate of the thresh- 
old amplitude T for the parametric instability to begin sap- 
ping energy from the eigenmode responsible for the standing 
wave resonance lock in our CO 6 simu lation from § |8| Us - 
ing the procedure detailed in § 6.5 of Bur kart et al.| ( 20T2] |, 
we found r~2x 10~ 8 . On the other hand, during the reso- 
nance lock the mode in question began with an amplitude of 
\q\ ~ 4 x 10~ 8 , which grew to ~ 10~ 7 before wave breaking 
destroyed the lock. 

This demonstrates that parametric instabilities may limit 
the achievable amplitudes of standing waves in close white 
dwarf binaries, potentially somewhat more stringently than 
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wave breaking alone. Exactly how this affects the overall tidal 
synchronization scenario requires more detailed study. 

10. CONCLUSION 

In this paper, we have studied the linear excitation of 
dynamical tides in white dwarf binaries inspiraling subject 
to gravitational wave radiation. We showed that the phe- 
nomenon of resonance locks occurs genetically in this sce- 
nario, both when the dynamical tide represents a standing 
wave or a traveling wave. ("Synchronization lock" may be 
a more appropriate term for this phenomenon in the trav- 
eling wave regime; however, we continue using "resonance 
lock" due to the many similarities between traveling wave and 
standing wave results.) 

In a resonance lock, as the orbital frequency increases ac- 
cording to ft = 57/fg W , where f gw is the gravitational wave in- 
spiral time (equation [TJ, a synchronizing torque produced by 
the dynamical tide causes the white dwarf spin frequency to 
evolve at nearly the same rate: f2 sp ; n rj f2 (§ kl. This means 
the l = m = 2 tidal driving frequency a = 2(fi^Tl sp j n ) remains 
constant, which in turn keeps the tidal torque nearly constant, 
leading to a stable situation. In other words, a resonance lock 
is a dynamical attractor (Figure[2]i. 

We first considered resonance locks created by standing 
waves, where resonances between the tidal driving frequency 
and white dwarf eigenmodes create the synchronizing torque 
required to maintain a ss constant. We derived analytic esti- 
mates of the orbital period P r \ at which such resonance locks 
can first occur (§|4j also Table [2]): P t \ ~ 30 min for hot car- 
bon/oxygen white dwarfs (T e jt ~ 12, 000 K) and P T \ ~ 200 min 
for cold carbon/oxygen white dwarfs (r e ff ~ 6,000 K). For 
helium white dwarfs, we found P r \ ~ 70 min for hot mod- 
els (T e H ~ 10,000 K), and P r i ~ 1 day for colder models 
(T e fi ~ 5,000 K). 

Tides preferentially deposit orbital angular momentum into 
a white dwarf's outermost layers, where wave damping is 
most efficient. A concern thus exists that a synchronously ro- 
tating critical layer might develop, causing rapid wave damp- 
ing and elimina ting the possibility of maintaining a standing 
wave (Goldreich & Nicholson 1989b). However, we showed 
that critical layers are in fact unlikely to develop in the stand- 
ing wave regime of white dwarf binary inspiral, since a typ- 
ical white dwarf fossil magnetic field is capable of winding 
up and enforcing solid-body rotation throughout t he w hite 
dwarf down to orbital periods of ~ 10 min or less (§ |6.2| Ap- 
pendix JA}. 

We derived analytic formulas for the tidal quality factor Q t 



(equation |l7[ ) and heating rate E^ e3t (equation 20 1 during a 
resonance lock (§ |5j. (Since Q t parametrizes the total tidal 
energy transfer rate, including mechanical energy associated 
with changing the white dwarf spin, values of Q t alone do 
not determine the tidal heating rate.) Our results predict that, 
for orbital periods of < hours, Q t ~ 10 7 for carbon/oxygen 
white dwarfs and Q t ~ 10 9 for helium white dwarfs. Our for- 
mula for Q t is independent of white dwarf eigenmode prop- 
erties and weakly dependent on the orbital period, scaling as 

—1 /3 

Qt oc f orb . It is, however, strongly dependent on the white 
dwarf mass and radius. We also found that tidal heating be- 
gins to rival typical white dwarf luminosities for P or b < 10 
min, a result that is relatively insensitive to white dwarf prop- 



population synthesis models for the evolution of close white 
dwarf binaries. 

As a standing wave resonance lock proceeds, the wave am- 
plitude required to maintain synchronization grows. Eventu- 
ally, the amplitude becomes so large that the stan ding wave 
begins to break near the surface convection zone (§ 6.1 1. This 



causes the dynamical tide to become a traveling wave, elimi- 
nating the resonance lock. This occurred soon after the initial 
resonance lock in our O.6M , T s s = 5,500 K carbon/oxygen 
white dwarf simulation; however, the standing wave reso- 
nance lock lasted much longer in our O.2M , r e f t = 7,000 K 
helium white dwarf simulation, from P or b ~ 250 min down to 
~ 40 min, amounting to ~ 10 Gyr of binary evolution. 

Resonance locks have traditionally been conside red only 
when the dy namical tide represents a standing wave ( Witte & 
|Savonije|199 9). We showed, however, that given sufficiently 
short orbital periods, reso nanc e locks can even occur in the 
traveling wave regime (§ 7.2 1. We derived two simple cri- 



teria for whether traveling waves can effect resonance locks: 
the traveling wave torque must be large enough to enforce 
f2 ss f2 sp j n , and the torque profile as a function of the tidal driv- 
ing frequency a = 2(f2-f2 sp j n ) must possess "jagged" features, 
a concept quantified by |iilogT/iilog(7| ^ 1 (equation 31 1, 
where r is the tidal torque. 

The first criterion is satisfied for orbital periods below a 
critical period P a \, which we found to be P a \ ^40-50 min in 
most white dwarf models (equation 33 Table [2}. The second 
criterion can be satisfied by rapid transitions between standing 
and traveling wave torques (which differ by orders of magni- 
tude) near resonances as a result of wave breaking (§ 6.1 1, or 
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The 



erties due to the steep power law scaling £heat 25 ^ c 
analytic results we derived can easily be incorporated into 



by wa ve in terference due to excitation by a composition gradi- 
ent (§ |7.1| Fulle r & Lai|2012| l. Excitation likely proceeds at a 
composition gradient in carbon/oxygen white dwarfs and hot 
helium white dwarfs, but excitation at the radiative-convective 
boundary becomes important for colde r hel ium white dwarfs 
with larger surface convection zones (§ |7.1| i. Excitation off a 
crystalline core may also be important in cold carbon/oxygen 
white dwarfs (§ |9.3| l. 

Even after the initial standing wave resonance lock is de- 
stroyed by wave breaking, a new traveling wave resonance 
lock takes hold once the orbital period declines to P or b ~ ^tri ~ 
40 - 50 min. The synchronization trajectory and correspond- 
ing values of the tidal quality factor (equation \Fi\ and tidal 
heating rate (equation [20| are the same during a traveling 
wave resonance lock. We confirmed our analytic derivations 
with numerical simulations that smoothly switched between 
standing and traveling wave torques based on the maximal 
value of the wave shear \k r ^ r \ (Figure B), wit h wa ve break- 
ing leading to traveling waves for | &,•<!;,• 1 (§ |6.1) . We pre- 
sented the results of two simulations, one with a helium white 
dwarf and one with a carbon/oxygen white dwarf, in §|8] Once 
the traveling wave resonance lock began, synchronization in 
our numerical calculations proceeded until the spin frequency 
f2 sp in became larger than 5fl = fi — S pin, meaning nonlinear 
rotational effec ts no t included in our analysis were likely to 
be important (§ |9.2| l. 

Our numerical calculations (Figure [6} demonstrate that ef- 
ficient tidal dissipation is produced by standing wave reso- 
nance locks at large orbital periods, and by traveling wave 
resonance locks at smaller orbital periods. The importance 
of the standing wave resonance lock regime at large orbital 
periods can be tested by measuring the rotation rates of wide 
white dwarf binaries. We predict that systems with orbital 
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periods of hours should have undergone significant synchro- 
nization (Figure [§}, while models that focus solely on exci- 
tation of traveling waves (Ful ler & Lai||2012] l would predict 
synchronization only at significantly shorter orbital periods. 
A second prediction of our model is that there may be a range 
of intermediate orbital periods (e.g., 20 min < P 01 -b < 40 min) 
where tidal dissipation is relatively inefficient compared to 
both smaller and somewhat larger orbital periods. 

The results derived here can be directly compared to the 
re centl y discovered 13-minute white dwarf binary J0651 
(§ 9.1 1. We predict a ~ 2% deviation of the orbital decay 
rate from the purely general relativistic value, which may be 
measurable given further observations. We also find that our 
predicted tidal heating rates are within an order of magnitude 
of the observed luminosities. This broad agreement is encour- 
aging given the well-known difficulties tidal theory has accu- 
rately predicting the efficiency of tidal dissipation in many 
stellar and planetary systems. 

In detail, we find that even if the helium white dwarf is 
nonrotating (which maximizes the tidal energy dissipated as 
heat), tidal heating is a factor of ~ 2 less than the observed lu- 



minosity, strongly suggesting that much of its luminosity must 
derive from residual nuclear burning or cooling of thermal en- 
ergy rather than tidal heating. In contrast, we predict that the 
carbon/oxygen white dwarf in J0651 should be ~ 10 times 
more luminous than is observed. We suspect that the origin 
of this discrepancy is the importance of rotational modifica- 
tion of stell ar ei genmodes at the short orbital period present 
in J0651 (§ 9.2 1, and perhaps the effects of nonlinear damp- 



ing/excitation of stellar oscillations (e.g. § 9.4 Weinberg et al. 
|2012| l. These will be studied in future work. 
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APPENDIX 



A. ANGULAR MOMENTUM TRANSPORT 

If we assume there is a source of angular momentum near the white dwarf surface, e.g. from tides, a fossil magnetic field of 
initial magnitude ~ Bo will wind up and exert magnetic tension forces attempting to enforce solid-body rotation. The rate at 
which magnetic tension transports polar angular momentum through a spherical surface S at radius r is given by 



1 

47T 



B r B& rsinddS, 



(Al) 



which can be derived by applying the divergence theorem to the magnetic tension force density (B • V )B/47r. As the field winds 
up, the radial component remains constant, while the azimuthal component increases as ( Spruit|1999| l 



Bd, =N vl B r , 



(A2) 



where /V w = rsin6* jidVL^/d^dt is the rotational displacement that occurs during wind up and we have assumed the rotational 

velocity field can be described by "shellular" rotation: v = ft sp i n (r)rsm9(j>. 

Once the field wind up propagates into the white dwarf core, an equilibrium field is established that communicates the tidal 
torque throughout the white dwarf and eliminates any rotational shear. There are several requirements necessary for this equilib- 
rium to be reached during inspiral, and consequently for white dwarfs in inspiraling binaries to rotate as a rigid bodies. 

A. 1 . Solid-body rotation at short orbital periods 

Whether there is sufficient time for the global equilibrium magnetic field to develop and eliminate differential rotation al- 
together amounts to whether the timescale over which the torque changes, given by the gravitational wave inspiral time f gw 
(equation fTJ, is longer than the timescale over which the wind up of the magnetic field propagates into the core, which is given 
by the global Alfven crossing time (?a) = / dr/v\, where va = B r / yfAitp is the radial Alfven speed. This restriction translates to 



where we have assumed B r ~ Bo ^ constant. We can evaluate this further as 
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The other requirement for solid-body rotation is that the rotational displacement /V w required for the equilibrium field config- 
uration mu st not be too extreme, since the field becomes susceptible to resistive dissipation as well as various instabilities as it 
winds up (Spru it|1999| . We can address this constraint by setting the tidal torque appropriate for a resonance lock r r = / r 17/f gw 



Linear tides in inspiraling white dwarf binaries 



17 



(equation |9j, where I r is the moment of inertia up to a radius r, equal to equation ( |A1[ >, and then requiring /V w < 1. Solving for 
Bo, we have 

I t o 

(A5) 

Since the right-hand side of equation (|A5|l scales radially as r 1 , we can safely set / f = and r = Evaluating this further yields 
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(A6) 



(A7) 



We see that equation ( |A4[ ) is likely to hold nearly until mass transfer, meaning there is always sufficient time to set up an 
equilibrium field capable of stably transmitting angular momentum throughout the white dwarf and enforcing solid-body rota- 
tion. Equation ( |A6[ > is more res trict ive, however, and shows that even for a fossil field of initial magnitude of Z?o ~ 10 6 G in 
a carbon/oxygen white dwarf (§ 6.2 1, excessive wind up may begin to occur below an orbital period of ~ 1 min. Nonethe less, 
at orbital periods of ~ 100 min where standing-wave resonance locks occur, which require solid-body rotation, equation ( |A6| l 
requires only modest fields of ~ 200 G. 

Lastly, although angular momentum transport occurs due to a wound-up equilibrium magnetic field, torsional Alfven waves 
can also be excited, which cause o scillations in the differential rotation profile. Formally, a phase-mixing timescale must elapse 
before such waves damp (Spruit 1999). However, due to the slow ly varying torque, we expec t Alfve n wave excitation to be 
weak. Indeed, in simulations of the solar magnetic field's evolution, Charbonneau & MacGregor (J 1993 1 found that Alfven wave 
amplitudes (in terms of <5f2 sp i n /J7 sp j n ) were small. 

A.2. Transport during an initial resonance lock 

When a standing wave resonance lock first begins to takes hold, angular momentum is applied exclusively to a thin layer near 
the outer wave turning point where wave damping predominantly occurs. Since no equilibrium field state has y et developed in 
this situation, a concern exists that the layer will rapidly synchronize and destroy the lock before it begins (§ |6.2) . 

To this end, we first compare the spin-up timescale of the layer to the Alfven travel time t& across it. Using the torque for a 
resonance lock from equation |9) and taking t& =H/v\, where H is the thickness of the layer, we have 



(2/3)4irpr 4 Hn s 

I* ^spin/^gw 



> 



B r 



Evaluating p at the radiative-convective boundary (RCB), letting B r ~ Bq, and setting rftiR, this becomes 



We can evaluate this further as 
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and p c is the central density. 

Secondly, even if no critical layer occurs, it is still necessary for the global, wound-up equilibrium field to develop quickly 
in order for an initial resonance lock to develop. Specifically, in order for a particular resonance to halt the increase in the tidal 
driving frequency a = 2(fl - fl s pkd> its locally applied tidal torque must be communicated globally fast enough relative to the 
timescale over whi ch th e torque changes significantly. If this cannot occur, the system sweeps through the resonance without 
locking (see also § |6.3| l. The timescale over which a maximally resonant torque decays by roughly a factor of two is given by 
the time for the detuning Su> to increase from zero to of order the associated mode's damping rate j n (see e.g. equation [5] also 
Appendix B.3 I. Comparing to the global Alfven travel time (t\), this condition becomes 



Bo> 




(All) 



Tn^gw JO 

which is identical to equation dA3b except with an additional factor of fi/7„ on the right-hand side. Evaluating further, we have 
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where Fj was defined in Appendix A. 1 The value of 100 G reduces to 6 G for a helium white dwarf (He7 from Table[lji, holding 
7„ constant. This i s a much stricter requirement than equation ( | A 1 0| >, but still seems likely to be satisfied given typical white 
dwarf fields (§E2l. 



B. GLOBAL NORMAL MODE ANALYSIS 
B . 1 . Mode dynamics 

Here we give an overview of linear normal mode analysis as it applies to tidal interactions. As such, we assume the fluid 
motions generated by the tidal potential represent standing waves; we discuss the possibility of traveling waves in §|7] 

In linear perturbation theory, we expand all fluid variables in spherical harmonics angularly, indexed by / and m, and adiabatic 
normal modes radially, indexed by the number of radial nodes n, as 

oo / 

SX(r, 0,c/>,t) = J2J2Yl qnim(t)5X nlm {r)Yi m (6,4>). (Bl) 

1=2 m=-l n 

We computed normal modes using the ADIPLS stellar pulsation package (Christensen-Dalsgaard 2008 ). In this work we concern 
ourselves only with quadrupolar eigenmodes, and consider only circula r orb its, meaning we can let / = 2 and m = ±2, since m = 
modes have no time dependence with zero eccentricity. Thus equation ( |B1[ ) becomes 

SX(rAM = 2^Re [q n (t)8X n (r)Y 12 {9 ^)] , (B2) 

n 

with m = 2 used throughout. 

We can write the momentum equation schematically as ( |Press & Teukolsky|1977[ ) 

p£ = - P £,[&-pVU, (B3) 

where the linear internal acceleration operator £ satisfies £[£„] = w„£ B , w„ is an eigenfrequency, and U is the tidal pote ntial. The 
corotating frame mode amplitude equations determining the behavior of each q„ c an be obtained from e quation ( |B3] > by taking 
the scalar product with £* on both sides and integrating over the star, yielding (e.g. Weinberg et al. 2012[ ) 

q„ + 2 ln q„ + U) 2 n q n = 2w 2 n eW Q n e iat , (B4) 

where W = ^3^/10, e = (M'/M)(R /af i s the tidal factor, a = m (Vt- f2 SD j n ) is the tidal driving frequenc y, an d the linear overlap 
integral Q„ is discussed in Appendix |B .4 ( Pres s & Teukolsky 1977 ). We have also inserted into equation ( |B4| i a damping rate 7„, 
the calculation of which we describe in Appendix B.3 We normalize our normal modes by setting 

E. =E„= f 2u 2 n (£„ + Kl + mln) P^dr. (B5) 
Jo 

Given slowly varying orbital and stellar properties, the steady-state solution to equation ( |B4| i is 

q„(t) = 2eQ„w(— 2 ^— )e-^'. (B6) 

The above amplitude applies in the corotating stellar frame; in the inertial frame the time dependence e~ ,a1 instead becomes 

B.2. Angular momentum and energy transfer 

Assuming a circular orbit and alignment of spin and orbital angular mo menta, th e secular quadrupolar tidal torque on a star is 
given by an expansion in quadrupolar (7 = 2) normal modes as (Burk artet al.|2012 Appendix CI and references therein) 

r = 8mE* eW £ Q \ 2 ^ 2 2 , (B7) 

where most variables are defined in the previous section. 
The tidal energy deposition rate into the star E^ t can be determined from r by the relation 

- Sorb = £tide = O T , (B8) 

valid only for a circular orbit. This can be derived by differentiating standard equations for the energy and angular momentum 
of a binary with respect to time, setting e = e = 0, and noting that the tidal torques and energy deposition rates in each star of a 

binary are independent. 

The total energy contained in the linear tide can be expressed in the corotating frame as ( Schenk et al.| 2002 ) 

Bade + ^ (B9) 
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where the operator £ was introduced in Appendix B. 1 Using results from Appendix B. 1 we can evaluate this expression as 



£,h = 2E e 2 W 2S To 2 " K 



(BIO) 



Lastly, since the great majority of the tidal energy is in the lowest-order modes, i.e. the equilibrium tide, which satisfy u n ^> a, 
we can further set a s» in the previous equation to derive the simple expression 



where 



£tide ~ A£*£ , 

\ = 2W 2 Y,Ql 



(Bll) 

(B12) 



Several limits can be taken of our general torque expression in equation ( |B7| i. First, if we assume to„ ^> (a, 7,,), we arrive at the 
equilibrium tide limit: 



T eq = 8mE,s 2 W 2 J2Q 2 n 



uj- 



(B13) 



|Willems et al.| ( |20"T0| showed that in linear theory the equilibrium tide provides a negligible torque due to the very weak damping 
present in white dwarfs, resulting from their high densities and long thermal times. We confirm this result: e.g., the tidal quality 
factor associated with damping of the quadrupolar equilibrium tide in white dwarfs is (equation 15 1 



Qeq _ \ Qn > i rA 

1 8£„e^7«K~ 



(B14) 



which is much larg er th an the effective value of Q t for the dynamical tide determined in § 5.1 The resonant dynamical tide 
(lo„ ~ a in equation B7 1 is discussed in §§|4]&[7] 

B.3. Damping 

We consider two different damping processes in this appendix: thermal diffusion and turbulent convection. We also calculate 
g-modes' group travel times, which sets the effective damping time for traveling waves. Figure [7] shows plots of all three of 
these quantities for several of our fiducial white dwarf models, and demonstrates that thermal diffusion is the dominant source of 
damping for high-order g-modes in white dwarfs. 

Damping due to electron conduction and radiative diffusion can be estimated simultaneously as ( |Goodman & Dickson|1998) 



in 



2.E n 



dr 
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where dE n /dr is the integrand of equation ( |B5[ ), \ is the thermal diffusivity (including both radiative diffusion and electron 
conduction), and k r is the radial wavenumber. In terms of an effective opacity n, \ can be expressed as 



I60T 3 

3np 2 c p 



For high-order g-modes, k r is given by ( Christensen-Dalsgaard 2008 
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where k\ = 1(1+ l)/r 2 is the angular wa venumber. The integration in equation ( |B15| l is performed up to the adiabatic cutoff radius 
defined by uit t h = 1 dUnno et al.|1989|, where f t h = pc p T/gF is the local thermal time. A wave's group travel time across a scale 
height must remain smaller than its local damping time in order for the wave to reflect. This criterion will always be broken for 
sufficiently long-period waves since the radial wavenumber k r grows as or 1 ; nonetheless, we find that thermal diffusion is never 
strong enough to invalidate the standing wave assumption for all of the modes used in this work that are capable of effecting 
standing wave resonance locks. 



To e stimate the turbulent convective damping rate, we rely on the calibration of convective viscosity performed by Penev et al. 
(2009). The formula we employ is 
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where sp' = 0.23 an d s\ = 0.084 (i Penev & Sasselov]|201 1 1. For the effective turbulent viscosity v im h, we use equation (11) of 
|Shiodeetal.| ( [20T2l ): 



f turb = Lv C onvmin 
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FIG. 7. — Plots of the contributions to eigenmode damping rates due to thermal diffusion 7* ff (thick blue lines) and turbulent convection 7'"* (red lines), as 
well as the effective damping rate for traveling waves a n = 2n /Jgroup,n (dashed green lines), as functions of the eigenmode frequency uj„ (in units of the dynamical 
frequency ) an d radial order n, for the first 400 g-modes in four of our fiducial white dwarf models (TablefTJ. We describe the computation of these quantities 
in Appendix |B.3| Inside the shaded regions in panels 2 & 4, g-modes become traveling waves due to rapid tnermal diffusion near their outer turning points. In 
this case, approximating wave damping rates using 7,f lff + 7, t , urb is invalid, and a traveling wave formalism such as that discussed in §^]must be used. Note that 
the white dwarf models shown here are different from those in Figure [4] 

where v conv is the convective velocity, L is the mixing length, ? ec idy = v conv /L, U m \ n = 0.1, and il max = 2.4. 

Lastly, the effective damping rate applicable in the traveling wave regime is the inverse group travel time a = 27r/f group ; we can 
calculate f gl - oup as 



^group — 2 



dr 



(B20) 



where v group = dtu jdk r and th e integration is over the propagation cavity where the wave frequency u> < N. Using the dispersion 
relation from equation ( |B17[ ), this becomes 
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For uj N, a oc ui 2 ; this proportionality is verified in Tableland Figure [7] 



B.4. Linear overlap integral 
The linear overlap integral for quadrupolar eigenmodes, introduced in equation dB4|, can be expressed as 
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with / = 2, where £h is the horizontal fluid displacement. The second equality can be derived by substituting the continuity 
equation, and the third equality by substituting Poisson's equation. However, all three of these methods of calculating Q„ suffer 
from numerical difficulties, presumably arising due to a failure of orthogonality o r completeness o f the numerically computed 
eigenmodes ( [Fuller & Lai|201 1) , or to small inconsistencies in the stellar model ( |Fuller & Lai|2012"] l. 

Thus we now consider a more stable way of numerically evaluating Q„, which is what we actually employed in our calculations 
and the fits in Table [3] We again focus on / = 2 modes, but the technique ca n ea sily be extended to arbitrary I. First, the tidally 
generated displacement field is give n as a sum of normal modes in eq uatio n ( |B1| >. If we set the tidal driving frequency a and the 
damping rate j n to zero in equation ( |B6[ > and substitute into equation ( |B2| i (while keeping the tidal factor e nonzero), we recover 
the equilibrium tide limit. However, the equilibrium tide can alternatively be obtained by directly solving the inhomogeneous 
linear stellar oscillation equations in the zero-frequency limit; see e.g. Weinberg et al7| ( |2012.| > Appendix A.l. Equating these two 
alternate expressions, taking the scalar product with p£* on both sides, integrating over the star, and solving for Q„, one obtains 

Q " = WkJ i X kC+Kl + ^^n)pr 2 dr, (B23) 

where 

X eq (r, e,<f>) = (X^(r)r+rX* (r) V) Y 2m (0, 0) (B24) 

m=±2 

is the numerically computed / = 2 equilibrium tide scaled to e = 1 . 

Finally, there is one further method of computing Q„ that we have employed, which also uses a solution to the inhomogeneous 
equations. In this method, however, instead of comparing alternate computations of the equilibrium tide, we instead evaluate the 
inhomogeneous tidal response very near an eigenfrequency. The overlap can then be extracted by fitting the resulting Lorentzian 
profile of e.g. the tidal energy. This method numerically agrees very well with the equilibrium tide method described above. 

C. WHITE DWARF MODELS 



We used MESA version 4298 ( Paxton et al.|201 1) to produce our helium white dwarf models (Table T| He5, He7, and Hel 0) 



We used three inlists. The first evolves a 1 .6M star with Z = 0.02 from ZAMS to where O.198M ofTts core is locally at least 
90% helium. Salient non-default parameter values are: 

mesh_delta_coef f = 0.5 
hl_boundary_limit = 0.1 
hl_boundary_mass_limit = 0.198 

The second smoothly removes the outer 1.4M Q , leaving a hydrogen layer with ~ 1% of the remaining mass. This is achieved 
with: 

relax_mass = .true. 
new_mass = 0.2 

The third evolves the resulting 0.2M Q white dwarf until T e ff = 5, 000 K. We invoke MES A's element diffusion routine even where 
the plasma interaction parameter T > 1, Salient non-default parameter choices are: 

mesh_delta_coef f = 0.05 

which_atm_option = ' simple_photosphere' 

surf_bc_of f set_f actor = 

do_element_dif fusion = .true. 

use_Ledoux_criterion = .true, 

dif fusion_gamma_full_of f = ldlO 

dif f usion_gamma_f ull_on = ldlO 

dif fusion_T_full_of f = -1 

dif f usion_T_f ull_on = -1 

dif fusion_Y_full_of f = -1 

dif f usion_Y_f ull_on = -1 

Our carbon/oxygen white dwarf models (Table [l] CO 6 and CO 12 ) were produced b y solving for hydrostatic equilibrium 
subject to heat transport by radiative diffusion and electron conduction (Hansen et al.|2004)>. We use the OPAL EOS and effective 



i (Ro^ 

opacities (Potekhin & Chabri er|2010[ ) where the electrons begin to become degenerate. We use a mixture of 25% carbon, 75% 
oxygen for the inner 98% of the model's mass, then add a helium layer with 1.7% of the mass, and finally a hydrogen layer 
with the remaining 0.17%. We smooth the composition transition regions over ~ Q.lH p with a Gaussian profile, where H p is a 



opacities (Rogers etal. 1996[) in the white dwarf outer layers, and transition to the Potekhin-Chabrier EOS and electron conduction 
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pressure scale height. We treat convection with mixing length theory, using L = H p for the mixing length. 
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